Full Counting Statistics for a Single-Electron Transistor, 
Non-equilibrium Effects at Intermediate Conductance 
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We evaluate the current distribution for a single-electron transistor with intermediate strength 
tunnel conductance. Using the Schwinger-Keldysh approach and the drone (Majorana) fermion 
representation we account for the renormalization of system parameters. Nonequilibrium effects 
induce a lifetime broadening of the charge-state levels, which suppress large current fluctuations. 
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The 'Full Counting Statistics' (FCS) of charge trans- 
port has proven to be a powerful tool in the description of 
current fluctuations 0. The concept had been explored 
by Levitov and Lesovik 0, who expressed the FCS of 
an arbitrary mesoscopic structure with non-interacting 
electrons in terms of its S— matrix. Much less is known 
about the FCS of interacting mesoscopic systems, a prob- 
lem which has been addressed only recently 0j0,0S0- 

As a fundamental example of an interacting meso- 
scopic systems we consider a Single-Electron Transistor 
(SET). It consists of a metallic island coupled to drain 
and source (left and right) electrodes via low-capacitance 
tunnel junctions, with resistances i?L and .Rr, as well as 
to a gate electrode. The strength of Coulomb interaction 
is characterized by the charging energy Ec = e 2 /2Cs, 
which depends on the total capacitance Cs = Cl+Cr+Cg- 
A measure for the tunneling strength is the dimcnsionlcss 
parameter ao = (i?L+-RR.)/27re 2 i?L-RR (we put H=kB = l)- 

In Refs. 0, 0| the FCS of a similar system - a quan- 
tum dot - has been studied, fully accounting for strong 
electron correlations, however only for a particular setup 
and parameters, corresponding to the Toulouse point. A 
renormalization group approach had been developed for 
the regime ao ^> 1 In the opposite limit, ao^O, the 
FCS has been analyzed to lowest order in Ref. [|| and 
next-to-lowest order (cotunneling) in Ref. 0. However, 
effects of quantum fluctuations induced by nonvanishing 
ao are still unknown. The aim of the present paper is to 
derive the FCS for a SET beyond perturbation theory in 
the intermediate strength tunneling regime ao<l. 

Let us further specify the situation to be considered. 
At low transport voltages and temperatures, eV, T<^Ec, 
due to Coulomb blockade tunneling is suppressed in a 
SET, everywhere except near specific values of the gate 
voltage, e.g., near Q G = CqVq = e/2. In the neighbor- 
hood of this conductance peak the Coulomb barrier is 
Ao = Ec(l — 2Qo/e). For a <C 1 electrons tunnel via 
the island sequentially only when iir < Ao < /xl, where 
A*l/r = k l/rzV is the voltage drop between the L/R elec- 
trode and the island, and kl/r = ±Cr/l(Cl + Cr)" 1 . 



With increasing ao, higher order effects such as cotun- 
neling and quantum fluctuations of the charge gain im- 
portance 0. They lead to a renormalization of Ao 
and an. The perturbative renormalization group anal- 
ysis yj (for eV — 0) predicts a renormalization factor 
zo = l/{l + 2ao \n(Ec / A) } to depend logarithmically on 
the cutoff energy A = max{Ao, T}. 

The model. - We concentrate on the tunneling regime 
with inverse RC time I/RtCs = ^na^Ec smaller than 
Ec, which ensures that the charge-state levels are well 
resolved. In the vicinity of the conductance peak, pre- 
cisely for Ao | /Eq -C 1, it is sufficient to restrict atten- 
tion to two charge states with charges differing by e. The 
Hamiltonian can then be mapped onto the 'multi-channel 
anisotropic Kondo model' |9|. Introducing a spin- 1/2 op- 
erator a acting on the charge states, we write 

H = Er=L,R,I J2kn £ rfc ® x kn® rk n. + ^t&z 

+ J2r=h,RJ2kk'n( T ^\ kn a r k'n^+ + H.C.). (1) 

Here a rk creates an electron with wave vector k and 
channel index (including spin) n in the left or right elec- 
trode or island (r=L,R,I). Tunneling matrix elements T r 
are assumed to be independent of k and n. The junction 
conductances are 1/R T = 27re 2 A c h|T r | 2 pip r , with 7V c h be- 
ing the number of channels and p r the electron DOS. We 
implicitly assume that energy and spin relaxation times 
are fast, and electrons obey the Fermi distribution. 

A convenient tool to treat the spin-1/2 operators in 
Eq. iJTJ is the 'drone' (Majorana) fermion representa- 
tion [lj, a + — c^4>, a z — 2&c — 1. where (f> = d) + d is a 
Majorana fermion and c and d are Dirac fermions. This 
formulation enables one to apply Wick's theorem and the 
fermionic Schwinger-Keldysh approach 0, 0] . 

Cumulant generating functional. - The central object 
of our approach is the generating functional of connected 
Green's functions (GFs) 



W[<p] = -i\nf D[a:. kn ,a rkn ,c*,c,d*,d}e^c^W. 



(2) 
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Here L is the Lagrangian corresponding to H QJ, and six 
Grassmann variables have been introduced (see Ref. 0] 
for details). The closed time-path C (Keldysh contour) 
runs from t = — oo to oo, back to — oo, and connects 
to the imaginary time-path to end at t — — oo — i/T. 
We introduce auxiliary source fields, the phase of the 
tunneling matrix element, T r — > T I e lKrip ( t \ distinguishing 
between forward and backward time-paths, <p+(t) and 
<P-(t). The 'center-of-mass' variable f c (t) = {<p+(t) + 
(y5_(t)}/2 = eVt is then fixed by the transport voltage. 

From the Cumulant Generating Functional (CGF) one 
finds the number of transmitted electrons q during the 
measurement time to, W(A) = X)^o((^9™))(*^)"/ n '- 
Following Ref. (lflj . it is derived from Eq. J3J by fix- 
ing during the measurement the 'counting field' (relative 
variable) (p&(t) = tp+(t) — tp—(t) at a constant value A: 

W(A) = iW[ip]\ iPcW=eVt ^ & (f )= xe(t /2+t)e(t /2-t)- (3) 

The distribution of q (or equivalently of the current I = 
eq/to) is given by the inverse Fourier transformation, 

p = — f d\e W(X) - lqX « e W(A*)-i(t //e)A^ / 4 n 

The integral can be evaluated in saddle point ap- 
proximation, with A* following from the relation / = 
— ied\W(X*)/t This approach is valid for long mea- 
surement times since W is proportional to to (see below) . 

We proceed following Ref. , where a conserving ap- 
proximation for the second cumulant had been developed. 
Tracing out the electron degrees of freedom leads to an 



effective action for the c and d fields, S = S c h 
composed of a charging and a tunneling term: 



S 



S ch = / dt{c(t)*(id t - A )c(t)+id(t)*d t d(t)}, (5) 
Jc 

S x = - [ dtdt 1 c*{t)<j)(t)a x (t,t')(j)(t')c(t') + 0{T?). (6) 
Jc 

Here a x — + a# is a particle-hole GF describing tun- 
neling of an electron from one electrode to the island. 
It depends on the counting field. The connection to the 
ordinary GF is established by a rotation by A r = n r X in 
the Keldysh space as follows: 



a x (uj) = U\ a T {u)U Xt , 5r(w) 



af(cu) 



where U\ T = exp(— iA r Ti/2) and [T%\ij = 1 — &y. The 
retarded and advanced components are given by 



{uj - /ir) 2 + E% ' 



where a r = l/2ire 2 R r , and the Keldysh component by 
af (w) = 2af (w)coth(cj - fi T )/2T. We introduced a 
Lorentzian cutoff to regularize the ultraviolet divergence 



and ignored the term 0(T" r 4 ) in the action ijBJ), since it is 
small in the limit ^> 1 . 

The free retarded GF of the Dirac fermion c, g** (u>) = 
l/(u) + if] — A ), describes the dynamics of charge exci- 
tations (77 is a positive infinitesimal). The corresponding 
self-energy E A = + accounts for quantum fluc- 
tuations of the island charge caused by tunneling. In- 
tegrating out G?-fields, we obtain the components of the 
self-energy in first order in olq: 



EfH=2af( W ), / — 



duo 1 ia«{uj') 



2ir uj + ir] — uj 1 



(8) 



(For simplicity we present here only the result for A = 0.) 
For a symmetric SET (R L = R R , C L = C R ), at T = and 
\uo\ < eV, one finds (u) » a \n(2E c /eV) w - ir/2, 
where T = Til + Tli + Tir + Tri is the sum of the rates 
r rI/Ir = ±(l/e 2 i? r ) (A -Mr)/(e ±(Ao " A1 ' )/T -l) describing 
tunneling into (out of) the island through the junction r, 
evaluated by Fermi's golden rule. 

We can proceed in a systematic diagrammatic expan- 
sion in ao 14]- To lowest order one obtains for the 
CGF: WW(A) = - J c dtdf g c (t,t')Y, x (t' ,t). We project 
the time from the Keldysh contour C to the real axis and 
observe that for long enough measurement times we can 

t /2 

approximate 5 to (w) = f_° t ' g , 2 dte~'' ut /2'7r by a J-function, 

S to (uj) -> 5(u), and (5 to {uj)) 2 -> t 5(u)/2TC. The latter 
ensures that any closed diagram, consequently W, is pro- 
portional to to- After Fourier transformation we obtain 

W [1] (A)w-t J dujTr{g c (Lu)T 1 Z x (Lo)T 1 }/2TT 

= t J2 {^- r ri(e iAr - 1) + P+T^e-^ - 1)}.(9) 

r=L,R 

Here we used the expression for the Keldysh component 
of c-field GF, g?(w) = 2i1mg?(w)(P--P+), which con- 
tains equilibrium occupation probabilities of the charge 
states Q G and Q G - e: P± = l/(e ±A °/ T + l). 

At this point we note that the naive first order ex- 
pansion @ is insufficient. First, it contains the equi- 
librium occupation probabilities rather than the station- 
ary ones. Second, due to charge conservation the CGF 
should depend only on the difference of the counting fields 
Al — Ar = A 0, which is also violated. These prob- 
lems are resolved if we sum up an infinite subclass of 
diagrams. Specifically, we sum up the geometric series in 
(<7 c TiE A Ti), which contains the leading logarithms, i.e. 
powers of a ln(2Ec/eV), and get 

W(A) =t [P Tr In k(w)- 1 - t 1 E^(o;)t 1 

/d 
— In [1 + T F {uj)f{uj - ^)h(u - mW X - 1) 



+ T F (uj)f(oj - nn)h(uj - (xt,)(e 



1)] 



(10) 
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Here /(w) = l/(e w / T + l), = 1 - /(«) and 

T» = -af ( W )af ( w )/| w - A - Sf H| 2 . (11) 

Note that we subtracted a constant from the CGF in 
order to satisfy the normalization condition W(0) = 0. 

Eq. (HHJ is the main result of this paper. It is similar 
to Levitov-Lesovik formula |2|, but the effective trans- 
mission probability (|11|) accounts for quantum fluctua- 
tions of the charge. Using the condition Eq. @, Eq. i|10|) 
can also be obtained from an approximate W given in 
Eq. (25) of Ref. Thus, the first and second cu- 

mulants, (7) = e((Sq))/t and S H = 2e 2 ((Sq 2 )) /t , 
reproduce the average current [lQ] and zero-frequency 
noise 0|, derived before. 

Although we used only the first order expansion for the 
self-energy, Eq. 1)11) [I is exact to second order in ao. One 
can check that the diagrams ignored in Eq. (|10|l within 
second order expansion, i.e. the diagrams with intersect- 
ing interaction lines, are proportional to Aq/Eq <C 1. 
Furthermore, higher order terms of Eq. (|10fl generate the 
renormalization factor zq consistent with the renormal- 
ization group result 0. 

Limiting cases. - In the limit ao — >0, Eq. 11U|I repro- 
duces the result of the 'orthodox' theory |(j : 



(12) 



d(x) = i + ^y-i) + 4r M y-i) (13) 

The second order expansion in ao reads, 

W (2) (A) = a Ao {RcSf (A )W (1) (A)} + W cot (A). (14) 

The first term of this equation provides the renormal- 
ization of the system parameters up to first order in 
ao. Namely, the parameters are renormalized as ao — > 
ao{l + <9 Ao Re£f (A )} and A -> A + Re£f (A ). This 
agrees with the corresponding results obtained earlier for 
the average current |l6(. It is also consistent with the 
recent results of Braggio et al. • We also checked that 
Eq. I|14|) can be reproduced by the systematic real-time 
diagrammatic expansion similar to that of Ref. . 

The second term of Eq. (|14|) is the CGF of a bidirec- 
tional Poissonian process 



W cot (A) =t {7 + (e lA -l)+7"(e 
governed by the cotunneling rates 



-iX 



1)}, (15) 



7 = / du) 



2ita};a$ (w - hl)(u - /i R ) 
(e± ^iL _ _ eT i^aj ™ ( W + lfl - Ao) 2 



Re- 



This term is relevant in the Coulomb blockade regime 
and is consistent with the FCS theory of quasiparticle 
tunneling [l7|. 
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FIG. 1: Current distribution P of a symmetric SET [T = 
and eV/Ec = 0.2]. (a) Solid lines are plots of P for various 
values of Ao. The dashed lines represent the 'orthodox' theory 
and the dot-dashed line represents the cotunneling expansion. 
The inset shows the average current for the same parameters. 

(b) Plot of P at Ao = for various values of the conductance 
as a function of the current normalized to V /2(JJl+JJr); inset: 
the same distribution normalized to the average current (7). 

(c) The skewness and (d) kurtosis for various conductances. 



At the conductance peak, Ao = 0, for a symmetric 
SET and T = 0, the 'orthodox' theory yields « 
2q{e lX / 2 - 1), where eq/t = V/2(R L + R n ). The fac- 
tor e lX > 2 leads to a sub-Poissonian value of the Fano 
factor Sn/2e(I) w 1/2, indicating that tunneling pro- 
cesses are correlated. The origin of this correlation can 
be understood from the explicit form of the distribution 
,q u =o Pp(qh)Pp(qR) S q ,( qij +q R )/2, obtained by 
inverse Fourier transformation of Eq. {2J without saddle 
point approximation. The numbers of electrons trans- 
mitted through either junction, <7l and <7r, follow the 
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same Poissonian distribution Pp(q) = q q e~^/q\. The 
Kronecker delta implies that and q^, are correlated. 

With A approaching the threshold value A /eV — 
0.5, the tunneling onto the island becomes the bottle- 
neck and the CGF acquires the Poissonian form, w 
to^hi (e lA — 1). It remains Poissonian in the cotunneling 
regime, |A /eV| > 0.5, W cot « i 7 + (e iX - 1). 

Let us compare our result (fT7Tf> to the 'orthodox' and 
the cotunneling theories. The latter theories fail around 
the threshold Ao/eV = 0.5. For example, in the aver- 
age current, shown in the inset of Fig. ^a), we observe 
a mismatch between their predictions since Tli — > 
while 7 + — > oo. In contrast, the distribution derived 
from Eq. shown in Fig. ^a) by solid lines, behaves 
regular. It widens with increasing Ao- The 'orthodox' 
(dashed lines) and the cotunneling (dot-dashed line) the- 
ories show the same trend but overestimate the width. 

Renormalization and lifetime broadening effects. - For 
large conductance, quantum fluctuations of the charge 
are pronounced. However, as long as zqT -C A, where 
A = max(|z A |,27rT, \eV\/2), the 'orthodox' CGF 
with renormalized parameters ZQCtQ and zqAq remains a 
good approximation. This scenario fails in the regime 
A <C Tk — Ece~ x l' laa /27r where the approximation of 
leading logarithms becomes insufficient. 

The renormalization effect is illustrated in Fig. ^b), 
where the current distribution at Ao = is plotted. 
Since Zo decreases with increasing ao, the mean value 
of the current, i.e. a peak position, shifts to lower values. 
The renormalization effect can be absorbed when we plot 
InP with horizontal axis normalized by (/) rather than 
y/2(i?L + Rr)- However, even after plotting the distri- 
bution as a function of the normalized current [inset of 
Fig. Gib)] the three curves do not collapse to a single 
one. The remaining differences can be attributed to the 
non-Markovian lifetime broadening effect as described by 
ImE^. We observe that the current distribution shrinks 
with increasing ctQ. This agrees with the previously noted 
suppression of the Fano factor [14j. FCS provides fur- 
ther information, showing in detail how the probability 
for currents exceeding the average value is suppressed. 

The effect of lifetime broadening is also visible in the 
moments. At moderately large voltages, eV 3> Tk, and 
at T = the real part of the self-energy is negligible 
and Sf (to) « -iira eV. The CGF at A = then is 

W(A) « 2g{(e lA / 2 -l)-2a (e lA -l) 

+ ir 2 a 2 (e l3X/2 - e lA/2 )/2 + O(ag)}, (16) 

and the ratio of higher order cumulants to the first one 
becomes ((Sq n )) / ((Sq)) = 2 1 ""{1 - 4a (2™- 1 - 1) + 
O(ag)}. We note that, as ao increases, higher order cu- 
mulants are suppressed as compared to the Poissonian 
result 2 1 -". 

Figures ^c) and (d) show the skewness ((5q 3 )) and 
the kurtosis ((Sq 4 )) as a function of A . A peak around 



the threshold develops with increasing conductance. We 
expect that this kind of behavior can be observed with 
present-day experimental techniques [l8j |. 

In conclusion, we have derived the Full Counting 
Statistics for a single-electron transistor in the vicinity of 
a conductance peak. Quantum fluctuations of the charge 
are taken into account by a summation of a certain sub- 
class of diagrams, which corresponds to the leading log- 
arithmic approximation. In first order in ao our results 
reproduce the 'orthodox' theory, while in second order 
they account for renormalization and cotunneling effects. 
We have shown that in non-equilibrium situations quan- 
tum fluctuations of the charge induce lifetime broadening 
for the charge states of the central island. An important 
consequence is the suppression of the probability for cur- 
rents larger than the average value. 
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